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Abstract 

The rational large eddy simulation (RLES) model is applied to turbulent channel 
flows. This approximate deconvolution model is based on a rational (subdiagonal Fade) 
approximation of the Fourier transform of the Gaussian filter and is proposed as an 
alternative to the gradient (also known as the nonlinear or tensor-diffusivity) model. 
We used a spectral element code to perform large eddy simulations of incompressible 
channel flows at Reynolds numbers based on the friction velocity and the channel half- 
width Rcr = 180 and Rcr = 395. We compared the RLES model with the gradient 
model. The RLES results showed a clear improvement over those corresponding to 
the gradient model, comparing well with the fine direct numerical simulation. For 
comparison, we also present results corresponding to a classical subgrid-scale eddy- 
viscosity model such as the standard Smagorinsky model. 
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1 Introduction 



Large eddy simulation (LES) is one of the most successful techniques in the numerical sim- 
ulation of turbulent flows. Contrary to the direct numerical simulation (DNS), which tries 
to capture all the scales in the flow, LES aims at resolving only the large-scale flow fea- 
tures. The large scales are defined by means of a filtering operation: the Navier-Stokes 
equations are convolved with a spatial filter, and the resulting filtered variables become the 
variables of interest in LES. Thus, a good LES model should be able to compute an accurate 
approximation of the filtered variables. 

An essential challenge in LES is the modeling of the subgrid-scale (SGS) stresses, rep- 
resenting the interactions between the large (above the filter width) and small (below the 
filter width) scales in the filtered Navier-Stokes equations. A remarkable research effort has 
led to a wide variety of SGS models, surveyed, for example, in 0, and 

Arguably the most popular class of LES models is the eddy-viscosity type, based on 
(variants of) the Smagorinsky model [Q. The main feature of the eddy- viscosity models is 
that they properly transfer kinetic energy (by inviscid processes) from large scales to smaller 
and smaller scales, until this energy is dissipated through viscous effects. These models have 
several limitations, however, including poor correlation coefficients in a priori tests |Q 
and inability to provide backscatter. Some of these limitations are circumvented by using a 
dynamic procedure in calculating the Smagorinsky constant, yielding the dynamic subgrid- 
scale eddy- viscosity model introduced by Germano et al. [0], and used in many studies [H, 

i- 

Another class of SGS models is the scale-similarity one. The scale-similarity model, 
introduced by Bardina et al. @], postulates that the full structure of the velocity field at 
scales below the filter width is similar to that at scales above the filter width. A priori tests 
show high correlations between real and modeled stresses. Another realistic feature of 
the scale-similarity model is that it produces backscatter. In a posteriori tests, however, 
the scale-similarity model does not dissipate enough energy and typically leads to inaccurate 
results. As a remedy, Bardina et al. added a dissipative Smagorinsky term. The resulting 
model, known as the mixed model, combines the strengths of both the scale-similarity and 
the Smagorinsky model. The dynamic procedure has been successfully applied to both the 
pure and the mixed scale-similarity model, yielding improved results [p!0| . 

A different class of SGS models consists of those models aimed at computing an improved 
SGS stress approximation by replacing the unknown unfiltered variables with approximately 
deconvolved filtered variables. An inverse filtered model was first proposed by Shah and 
Ferziger [|1^]. This idea was formalized by Geurts |T^ for the top hat filter. Kuerten et 



al. [0 used the approximate inverse to improve the computable estimates in the dynamic 
Smagorinsky model. Another model in this class is the velocity estimation model of Do- 
maradzki and Saiki [|^, [|15[, [1^. Stolz and Adams fl^ developed the approximate decon- 



volution model, based on repeated application of the filter to approximately deconvolve the 
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dependent variables [13], [p 



One popular model in this class is the gradient model (also known as the nonlinear or 
tensor-diffusivity model), which uses explicit filtering. In addition to the implicit filtering 
due to the effective truncations (grid and numerical method), this LES model also assumes 
a regular explicit filter of prescribed shape and effective width larger than the grid spacing. 

The gradient model is based on a Taylor series approximation of the Fourier transform 
of the filter and aims at reconstructing the filtered-scale stress due to explicit filtering. The 
gradient model was developed in several steps. First, in 1974 Leonard pO[ proposed a model 
for the "resolved scales" u u in the Reynolds stress tensor. Next, in 1979 Clark, Ferziger, 
and Reynolds used the same approach to model the "cross-terms" uu' + u'u. 

The gradient model was tested a priori against experimental data (two-dimensional cuts) 
by Liu et al. Borue and Orszag |^ presented a detailed a priori analysis of the gradient 
model based on Gaussian-filtered DNS of homogeneous, isotropic decaying turbulence. Also, 
Winckelmans et al. |]23[ presented several a priori tests for the gradient model and its dynamic 
version, again in the context of homogeneous, isotropic decaying turbulence. Similar tests 
have been performed by Carati et al. |^. All the above a priori tests have shown high 
correlations. 

In a posteriori tests, however, it was found that the gradient model does not dissipate 
enough energy. Simulations with the pure gradient model appear to be unstable [^. Also, 
Liu, Meneveau, and Katz |^T| reported problems near the wall, where the pure gradient 
model's Reynolds stresses do not follow the behavior. To stabilize the gradient model, 
Clark, Ferziger, and Reynolds combined it with a Smagorinsky term, but the resulting 
mixed model inherited the excessive dissipation of the Smagorinsky model. A different 
approach was proposed by Liu et al. 
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who supplied the gradient model with a "limiter" 
to prevent energy backscatter; this clipping procedure ensures that the model dissipates 
energy from large to small scales. This approach was also used in |2^, p7| . 

From this point of view, the gradient model is similar to the scale-similarity model: it 
shows high correlations in a priori tests, but it does not dissipate enough energy in actual 
LES simulations: hence the need for extra viscosity type terms (mixed models). We note 
that, for both types of model, the best results in actual LES simulations were obtained by 

I, i, m, 



using the dynamic mixed procedure |28|, |23 



In fact, it has been noted before [pq] , 
that there are strong ties between the gradient model and the scale-similarity model: the 
first term in the Taylor series expansion of the scale-similarity model is indeed the gradient 
model. As noted by Winckelmans et al. |23|, however, the other terms in the expansion are 



different. Thus, the gradient model is not identical to the scale-similarity model. 



The model presented in this paper was introduced by Galdi and Layton p9[ as an al- 
ternative to the gradient model. They observed that the Taylor series approximation of the 
Fourier transform of the Gaussian filter used in the derivation of the gradient model actually 
increases the high wave number components, instead of damping them. As an alternative to 
the Taylor series approximation, Galdi and Layton proposed a rational ((0,1) Fade) approx- 
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imation. This rational approximation is consistent with the original approximated function 
(which is a negative exponential): it attenuates the high wave number components. 

In this paper, the resulting LES model, called in the sequel the rational LES (RLES) 
model, is applied to the numerical simulations of incompressible channel flows at Rcr = 180 
and Rcr = 395. 



2 The Rational LES Model 

The usual LES starts by convolving the Navier-Stokes equations (NSEs) with a spatial filter 
gs- Assuming that differentiation and convolution commute (which is true for homogeneous 
filters), the filtered NSEs read as follows: 

+ V • (uu) - i?e"^Au + Vp = f , (1) 

where 5 is the filter width and u = * u is the variable of interest. The filtered NSEs (0) 
do not form a closed system, and a considerable research effort in LES research has been 
directed at modeling the stress 

r = uu — u u. (2) 

As mentioned by Carati et al. [|^, this stress consists of a filtered-scale stress tensor, 
mainly due to filtering, and a subgrid-scale (SGS) stress tensor, mainly due to discretization. 
One way of approximating the filtered-scale stress tensor is by using a Taylor series expansion 
in the wave number space to represent the unknown full velocity in terms of the filtered 
velocity. This approach was first used by Leonard [^, and it was later espoused by Clark, 
Ferziger, and Reynolds f^. The resulting model, called the gradient, nonlinear, or tensor- 
diffusivity model, was used in numerous studies g,§,[|g,|2|, [g,|g,|2T|,|T 



The gradient model is derived by using a Taylor series approximation to the Fourier 
transform of the Gaussian filter 

gsik) = ^ 1 - ^ + 0(5^), (3) 



and for its inverse 
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e— ^ l + lli^ + 0(5^). (4) 



^5(k) 47 
Decomposing u into its average and its turbulent fluctuations 

u = u + u', (5) 
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and taking first the average and then the Fourier transform of the above relation, we get 



S^(^-i)n. (6) 

and thus, 

u = 4u, (7) 

95 

where u denotes the Fourier transform of u. 

By taking the inverse Fourier transform and using (^), we get 

u ^ u H An. (8) 

47 ^ ' 

By plugging the above into (0), using (|^) and the same technique as above, simplifying, and 
dropping out the terms of 0(5^), we get the gradient model 

__ 52 

r = uu — u u ^ — VuVu, (9) 
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where 

d 

dxi d^i 



(VuVuK, = ^^^. (10) 

1=1 



Noticing that the approximation by Taylor series of gs actually increases the high wave 
number components (see Figure |l|), Galdi and Layton developed a new LES model based 
on a rational ((0,1) Fade) approximation of gs, which preserves the decay of the high wave 
number components: 

^^(k) = e-'-^ ^ + 



1 + 



47 

The resulting LES model, called the rational LES (RLES) model, reads as follows: 

-|:A + /)"(|1vuVu)]. (12) 

The inverse operator in (|12D acts as a smoothing operator and represents the approxima- 
tion of the convolution by the Gaussian filter in the stress tensor r in (□). 
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We note that differential filters have been proposed by Germano in |Q: Actually, one 
can think of ([T^) as the stress tensor obtained by applying such a differential filter. Mullen 
and Fischer used similar filters in [Bl]. Also, Domaradzki and Holm considered the Navier- 
Stokes-alpha model (which contains an inverse operator similar to the one in (p!2D), in an 
LES framework [p5[] . 

The mathematical analysis associated with the RLES model (p!2D was presented in [p6| ; 
the smoothing property of the inverse operator in (|T^) eliminated the necessity for using 



additional regularization operators (of eddy- viscosity type), as for the gradient model [p7 
The first steps in the numerical analysis and validation of the RLES model (|T 
and 1^ 



were made 



m 
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respectively. 

This paper presents numerical results for the RLES model (|l^) applied to the 3D channel 
flow test problem at Reynolds numbers based on the wall shear velocity Rcr = 180 and 
Rcr = 395. Some preliminary work started in it was significantly updated and improved 
in the present paper. 



3 Numerical Setting 

The 3D channel flow (Figure 0) is one of the most popular test problems for the investigation 



of wall bounded turbulent flows |^T|, [^]. We used the fine DNS of Moser, Kim, and Mansour 
P3| as benchmark for our LES simulations. 
We compared the RLES model (|12]) with 

(I) the gradient model (H) r = |;^VuVu; 

(II) the Smagorinsky model r = — ((7^5)^ l^l S; 

(III) a coarse DNS (no LES model), 

where S := |(Vu + Vu"^) is the deformation tensor of the filtered field, Cg = 0.1 is the 
Smagorinsky constant, and 7 = 6 is the parameter in the definition of the Gaussian filter. 

The computational domain is periodic in the streamwise (x) and spanwise (z) directions, 
and the pressure gradient that drives the flow is adjusted dynamically to maintain a constant 
mass flux through the channel. The parameters used in the numerical simulations are given 
in Table ^for the two Reynolds numbers considered {Rcr = 180 and Rcr = 395). The filter 
width 6 is computed as 6 = ^ Ax Ay{y), where A^^ and A^ are the largest spaces between 
the Gauss-Lobatto-Legendre (GLL) points in the x and z directions, respectively, and Ay{y) 
is inhomogeneous and is computed as an interpolation function that is zero at the wall and 
is twice the normal mesh size for the elements in the center of the channel. 

The numerical simulations were performed by using a spectral element code based on 
the Ptv — P N-2 velocity and pressure spaces introduced by Maday and Patera |Q . The 



domain is decomposed into spectral elements, as shown in Figure Mesh spacing in the 
wall-normal direction (y) was chosen to be roughly equivalent to a Chebychev distribution 
having the same number of points. The velocity is continuous across element interfaces and 
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is represented by A^th-order tensor-product Lagrange polynomials based on the GLL points. 
The pressure is discontinuous and is represented by tensor-product polynomials of degree 
N — 2. Time stepping is based on an operator-splitting of the discrete system, which leads to 
separate convective, viscous, and pressure subproblems without the need for ad hoc pressure 
boundary conditions. A filter, which removes 2%-5% of the highest velocity mode, is used to 



stabilize the Galerkin formulation |45]; the filter does not compromise the spectral accuracy. 



Details of the discretization and solution algorithm are given in [^], |^ . 

The initial conditions for the Rcr = 180 simulations were obtained by superimposing a 2D 
ToUmien-Schhchting (TS) mode of 2% amplitude and a 3D TS mode of 1% amplitude on a 
parabolic mean flow (Poiseuille flow) and integrating the flow for a long time (approximately 
200 B./ur, where H is the channel's halfwidth and Ur is the wall shear velocity) on a finer 
mesh (72 x 73 x 72 mesh points). The final field file was further integrated on the actual 
coarse LES mesh (36 x 37 x 36 mesh points) for approximately 50 H/m^ to obtain the initial 
condition for all four Rcr = 180 simulations. 

The initial condition for the Rcr = 395 case was obtained in a similar manner: We 
started with a field file corresponding to a Rcr = 180 simulation, and we integrated it on a 
finer mesh (96 x 73 x 72 mesh points) for a long time (approximately 50 H/u^). Then, we 
integrated the resulting flow on the actual coarser LES mesh (72 x 55 x 54 mesh points) for 
another 40 B./ur, and the final field file was used as initial condition for all four simulations. 

For each of the four simulations and for both Rcr = 180 and Rcr = 395, the flow was 
integrated further in time until the statistically steady state was reached (for approximately 
15 H/ttr). The statistically steady state was identified by a linear total shear stress profile 
(see Figure § and Figure H). 

The statistics were then collected over another 5 H/m^ and contained samples taken after 
each time step {At = 0.0002 for Re^ = 180 and At = 0.00025 for Re^ = 395). We also 
averaged over the two halves of the channel. 

In our numerical experiments, we considered, as a first step, homogeneous boundary 
conditions for all LES models tested. 

The numerical results include plots of the following time- and plane-averaged quanti- 
ties normalized by the computed Ur- the mean streamwise velocity <^ u ^ /ur] the x^y- 
component of the Reynolds stress <^ u'v' ^ /m^; and the rms values of the streamwise 
<^ u'u' ^ /u1; wall-normal -C v'v' ^ /m^, and spanwise -C w'w' S> /m^ velocity fluctua- 
tions, where <^ ■ ^ denotes time and plane (xz) averaging, the fluctuating quantities /' are 
calculated as /' = /— ^ / and a superscript denotes the variable in wall units. 

Note that in our simulations the bulk velocity Um was fixed to match the corresponding 
one in |4^ (see Table and the friction velocity u-r was a result of the simulations. Table 2 
presents the actual values of Rcr corresponding to the friction velocity Ur computed for all 
four numerical tests and two nominal Reynolds numbers. We note that the friction velocity 
Ut is within l%-2% of the nominal value, and, as a result, so is the actual Rcr- 
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4 A Posteriori Tests for Rcr = 180 



We ran a posteriori tests for the RLES model (12), the gradient model @, the Smagorinsky 



model, and a coarse DNS (no model). We compared the corresponding results with the fine 
DNS simulation of Moser, Kim, and Mansour |^ 



Figure ^ shows the normalized mean streamwise velocity; note the almost perfect over- 
lapping of the results corresponding to the models tested. We interpret this behavior as a 
measure of our success in enforcing a constant mass fiux through the channel. 

Figure |^ presents the normalized x, y-component of the Reynolds stress. The RLES 
model is a clear improvement over the rest (i.e., closest to the fine DNS). 

Similarly, the normalized rms values of the streamwise velocity fiuctuations in Figure 
^ show a better (closer to the fine DNS benchmark results in |^) behavior for the RLES 
model dH). 



Figures y and |T0|, containing the normalized rms values for the wall-normal and spanwise 
velocity fiuctuations, merit a more detailed discussion. Figure ^ shows the normalized rms 
values of the wall-normal velocity fiuctuations. Here, the RLES model (p!2| ) performs worse 
than the gradient model (P) near the wall. Away from the wall, the gradient and the RLES 
model perform similarly. The best results are obtained with the Smagorinsky model. 

The normalized rms values of the spanwise velocity fiuctuations in Figure ^ are better 
for the RLES model (|T2p than for the gradient model, except for a portion of 1 — \y\ roughly 
between 0.3 and 0.5. The Smagorinsky model gives the best results near the wall, but it 
overpredicts the correct value near the center of the channel. 

In conclusion, the RLES model ([T^ ) performs better than the gradient model, with the 
exception of the normalized rms values of the wall-normal velocity fiuctuations. The RLES 
model (jT^) is also more stable numerically than the gradient model. 



5 A Posteriori Tests for Rcr = 395 



We ran simulations with all four models for Rcr = 395, and we compared our results with 



the fine DNS in 03]. Again, as in the Rcr = 180 case, the normalized mean streamwise 



velocity fluctuations in Figure ^ are practically identical; this time, however, they do not 
overlap onto that corresponding to the fine DNS. Nevertheless, the mean fiows are the same, 
and this is supported by the fact that the models underpredict the correct value near the 
wall but overpredict it away from the wall. 



The normalized x, y-component of the Reynolds stress in Figure 12 is almost identical 



for all four models. This behavior was also noticed by Winckelmans et al. ||23[| . 

The same behavior can be noticed for the normalized rms values of the streamwise velocity 
fiuctuations in Figure |l^: the profiles for the four models are almost identical. The RLES 
model (|T^) performs slightly better near the wall, and the gradient model performs slightly 
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better away from the wall. They both perform better than the Smagorinsky model near the 
center of the channel. 



Figure [T3 presents the normalized rms values of the wall-normal velocity fluctuations. 
As for the i?e^ = 180 case, the gradient model performs better near the wall, and the RLES 
model (|12D performs better away from the wall. The Smagorinsky model performs best, but 
it overpredicts the correct value near the center of the channel. 

The same behavior is observed for the normalized rms values of the spanwise velocity 
fluctuations in Figure and again the Smagorinsky model performs best. 

In conclusion, for the Rcr = 395 case, the gradient and the RLES model ([121) yield 
comparable results. The best results, however, are obtained by using the Smagorinsky 
model. 

Again, as in the Rcr = 180 case, the RLES model ( ]T^ ) is much more stable numerically 
than the gradient model. 



6 Conclusions 

We have used a spectral element code to test the RLES model (|T^) in the numerical sim- 
ulation of incompressible channel flows at Rcr = 180 and Rcr = 395. This approximate 
deconvolution model is based on a rational (Fade) approximation to the Fourier transform 
of the Gaussian fllter and is proposed as an alternative to the gradient model (|^). We com- 
pared the RLES model ([T2|) with the gradient model, the Smagorinsky model, and a coarse 
DNS (no LES model). The corresponding results were benchmarked against the flne DNS 
calculation of Moser, Kim, and Mansour |H3l . 



The RLES model (p!2D yielded the best results for the Re^ = 180 case. These improved 
results were accompanied by a much increased numerical stability compared with the gradient 
model. 

The situation was different for the Rcr = 395 case. Here the RLES model ( p!^ ) and the 
gradient model yielded comparable results, and the Smagorinsky model performed the best. 
Again, the RLES model (|T2p was much more stable numerically than the gradient model. 

We beheve that these results for the RLES model are encouraging. They also support our 
initial thoughts: The RLES model is an improvement over the gradient model as a subfllter- 
scale model, and this is illustrated by the improved results for the Rcr = 180 case. The 
RLES model is also more stable numerically because of the additional smoothing operator, 
and this feature is manifest for both low [Rcr = 180) and moderate {Rcr = 395) Reynolds 
number flows. 

However, the RLES model accounts just for the subfllter-scale part of the stress recon- 
struction. The information lost at the subgrid-scale level must be accounted for in a different 
way, as advocated by Carati et al. |2^. This was illustrated by the fact that, even for a mod- 
erate Reynolds number {Rcr = 395) flow, the Smagorinsky model, a classical eddy-viscosity 
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model, performed best. 

Along these lines, our next step will be to develop a mixed model, consisting of the RLES 
model supplemented by a Smagorinsky model. We also plan to study improved boundary 
conditions, the commutation error [Q, [^], and the relationship between the filter radius 



and the mesh-size in a spectral element discretization. 
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Figure 1: Approximations to the Fourier transform of the Gaussian filter: Rational (Fade) 
vs. Taylor. 
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Table 1: Parameters for the numerical simulations. 



Nominal iie^ 


X Ly X L^: 




X 


Ny 


X 




180 


47r X 2 X Itt 


36 


X 


37 


X 


36 


395 


27r X 2 X TT 


72 


X 


55 


X 


54 
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Figure 2: Problem setup for the channel flow. 
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Figure 3: Spectral element meshes: Rcr = 180 (top), and Rct = 395 (bottom). 
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Table 2: Computed Ur and Rcr. 



Fixed Um Nominal Re-r 



Case 



Computed Ur Computed Re-r 



15.63 



180 



RLES 
gradient 
Smagorinsky 
coarse DNS 



0.9879448 
0.9890118 
0.9917144 
0.9873800 



177.8352 
178.0222 
178.5120 
177.7184 



17.54 395 RLES 1.001025319 395.4071960 

gradient 1.005021334 396.9859924 

Smagorinsky 0.9974176884 393.9718933 

coarse DNS 0.9901855588 391.1294861 



21 




Figure 6: Mean streamwise velocity, Rcr = 180. We compared the RLES model ([I2D, the 
gradient model (]^), the Smagorinsky model, and a coarse DNS, with the fine DNS of Moser, 
Kim, and Mansour P3 |. 
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Figure 7: The x, y-component of the Reynolds stress, Re^ = 180. We compared the RLES 
model (|T^, the gradient model the Smagorinsky model, and a coarse DNS, with the 
fine DNS of Moser, Kim, and Mansour ||43|| . 
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Moser et al. 




1-|y| 

Figure 8: Rms values of streamwise velocity fluctuations, Re-r = 180. We compared the 
RLES model (0), the gradient model the Smagorinsky model, and a coarse DNS, with 
the flne DNS of Moser, Kim, and Mansour [Q. 
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Figure 9: Rms values of wall-normal velocity fluctuations, Rer = 180. We compared the 
RLES model (0), the gradient model the Smagorinsky model, and a coarse DNS, with 
the fine DNS of Moser, Kim, and Mansour [Q. 
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Figure 10: Rms values of spanwise velocity fluctuations, Re-r = 180. We compared the 
RLES model (0), the gradient model the Smagorinsky model, and a coarse DNS, with 
the fine DNS of Moser, Kim, and Mansour [Q. 
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Figure 11: Mean streamwise velocity, = 395. We compared the RLES model ([T^), the 
gradient model (]^), the Smagorinsky model, and a coarse DNS, with the fine DNS of Moser, 
Kim, and Mansour P3 |. 
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Figure 12: The x, y-component of the Reynolds stress, Rcr = 395. We compared the RLES 
model (|T^, the gradient model the Smagorinsky model, and a coarse DNS, with the 



fine DNS of Moser, Kim, and Mansour |^ 



28 




Figure 13: Rms values of streamwise velocity fluctuations, Rcr = 395. We compared the 
RLES model (0), the gradient model the Smagorinsky model, and a coarse DNS, with 
the fine DNS of Moser, Kim, and Mansour 
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i-|y| 

Figure 14: Rms values of wall-normal velocity fluctuations, Rcr = 395. We compared the 
RLES model (0), the gradient model the Smagorinsky model, and a coarse DNS, with 
the fine DNS of Moser, Kim, and Mansour [Q. 
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Figure 15: Rms values of spanwise velocity fluctuations, Re-r = 395. We compared the 
RLES model (0), the gradient model the Smagorinsky model, and a coarse DNS, with 
the fine DNS of Moser, Kim, and Mansour [Q. 
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